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0\ . Abstract 



In this paper we employ recent results from real algebraic geometry and theory of 
moment problems to make the first step towards resolving the question of existence of 
I multivariate tight wavelet frames whose generators have at least one vanishing moment. 

The so-called Unitary Extension Principle from [M] and the results in [32] allow us to re- 
formulate the question of existence of tight wavelet frame in terms of the existence of the 
sum of squares decomposition of a single trigonometric polynomial with real coefficients. 
Our main result confirms the existence of such decompositions in the two-dimensional 
case. We also give sufficient conditions for existence of tight wavelet frames in the di- 
mension d >3 and illustrate our results with several examples. 
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1 Introduction 

(N 

There are several fundamental results by two groups of authors I. Daubechies, B. Han, A. Ron, 
Z. Shen [17] and C. Chui, W. He, J. Stockier |9l HD] that build up the theory of tight wavelet 
^ . frames and also provide their characterizations. Those characterizations, on the one hand, al- 
^ ! low us to establish the connection between frame constructions and a difficult algebraic problem 
of existence of sums of squares representations (sos) of non-negative trigonometric polynomials. 
On the other hand, these characterizations also offer methods, however unsatisfactory from 
the practical point of view, for construction of tight wavelet frames. 

The existence and practical methods for construction of tight frames, together with good 
estimates on the number of frame generators, are still open problems. On the one hand, one 
could get discouraged by a general result by Scheiderer in [37], which implies that not all 
non-negative trigonometric polynomials in the dimension d > 3 possess sos representations. 
On the other hand, we are dealing mostly with the case d = 2 and non-negative trigonometric 
polynomials with special properties. This motivates us to pursue the issue of existence of sos 
representations further. 
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It has been observed in [15] that redundancy of wavelet frames has advantages for apphca- 
tions in signal denoising - if the data is redundant, then loosing some data during transmission 
does not necessarily effect the reconstruction of the original signal. Shen et al. ^Jlj use the 
tight wavelet frame decomposition to recover a clear image from a single motion-blurred image. 
In [T] the authors show how to use multiresolution wavelet filters P and Qj to construct ir- 
reducible representations for the Cuntz algebra and, conversely, how to recover wavelet filters 
from these representations. Wavelet and frame decompositions for subdivision surfaces are 
one of the basic tools, e.g., for progressive compression of 3-d meshes or interactive surface 
viewing [5l EHl [33]. Adaptive numerical methods based on wavelet frame discretizations have 
yielded very promising results [I2l [13] when applied to a large class of operator equations, in 
particular, PDE and integral equations. 

In this paper we employ recent results from real algebraic geometry and theory of moment 
problems to make the first step towards resolving the question of existence of multivariate 
tight wavelet frames with generators having at least one vanishing moment. In other words 
the foundation for our results is the so-called Unitary Extension Principle from [51], a special 
case of the above mentioned characterizations in [ni[ini[IZ]- In section [31 we give several equiv- 
alent formulations of UEP that allow us to reformulate the problem of construction of tight 
wavelet frames as a problem of semi-definite programming. This establishes a connection be- 
tween constructions of tight wavelet frames and moment problems, see [2H |29l [30] for details. 
In section [JTl using results of [SH], we show that the existence of such tight wavelet frames 
in the two-dimensional case is always guaranteed. We also present an elegant method that 
sometimes simplifies the actual frame construction and illustrate this method on the example 
of the so-called butterfly scheme from [20] . In section \A.2\ we give a sufficient condition for 
the existence of tight wavelet frames in dimension d > 3 and illustrate our results with several 
examples of three-dimensional subdivision. We also show that not all 3-dimensional trigono- 
metric polynomials under consideration can be written as sums of squares of trigonometric 
polynomials with real coefficients. 

We list some existing constructions of compactly supported MRA wavelet tight frames 
of L2(M°') [H [m im [231 [321 [31 [38] that employ the Unitary Extension Principle. For any 
dimension and in the case of a general expansive dilation matrix, the existence of tight wavelet 
frames is always ensured by [21 [1], if the coefficients of the associated refinement equation are 
real and nonnegative. There are only few compactly supported multi-wavelet tight frames in 
the hterature, see [3l[5l[22]. In [3] the authors present a general method for constructing multi- 
wavelet tight frames in the case when the matrix coefficients of P all have only non-negative 
entries. The results of [S] are based on the time-domain techniques developed in [51 [ID] and boil 
down to symmetric factorizations of local positive semi-definite real matrices. The method in 
[5] for matrix-valued P is applicable only in the special case, when the refinement coefficients 
satisfy identity [5, (9)]. Another construction approach is given in ^2] and is an adaptation of 
the method based on singular-value decompositions first presented in [31] . 
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2 Background and Notation 



Let d e N, let T denote the d-dimensional anisotropic real (algebraic) torus, and let M[T] 
denote the (real) affine coordinate ring of T 

R[T] = R[xj,yf.j = l,...,d]/{x^j+y]-l:j = l,...,d). 

Rather than working with this description, we will mostly employ the complexification of T, 
together with its affine coordinate ring C[T] = M.[T] (8)r C. This coordinate ring comes with a 
natural C/R- involution * on C[r], induced by complex conjugation. Namely, 

C[T] = C[zt\...,zr] 

is the ring of complex Laurent polynomials, and * sends Zj to zj^ and is complex conjugation 
on coefficients. The real coordinate ring ]R[T] consists of the *-invariant polynomials in C[T]. 

The group of C-points of T is T(C) = (C*)'^ = C* x • • • x C*. In this paper we often denote 
the group of R-points of T by T*^. Therefore, 

T-^ = T(M) = {{z,,...,za)e{C*)'': = ■ ■ ■ = |^,| = 1} 

is the direct product of d copies of the circle group S^. The neutral element of this group we 
denote by 1 = (1, ... , 1). 

Via the exponential map exp, the coordinate ring C[T] = C[zf,..., 4] of T is identified 
with the algebra of (complex) trigonometric polynomials. Namely, exp identifies Zj with e~^'^^ , 
j — 1, . . . ,d. In the same way, the real coordinate ring R[T] is identified with the ring of 
real trigonometric polynomials, i.e. polynomials with real coefficients in cos(a;j) and sm{ujj), 
j = l,...,d. 

Let M e Z"'^'^ be a general expansive matrix, i.e. p{M~^) < 1, or equivalently, all eigen- 
values of M are strictly larger than 1 in modulus. Let m := |det M\. The finite abelian group 
G := 2ttM-'^Z'^/2ttZ'^ has order |G| = m, its character group is G' = Z'^/M'L'^. 

A wavelet tight frame is a structured system of functions that has some special group 
structure and is defined by the actions of translates and dilates on a finite set of functions 
V'j £ -^2(R^), ^ ^ j ^ N . More precisely, we define translation operators on L2(M"') by 
T„/ = /(■ - a), a G Z'', and dilation (homothethy) Um by UmI = m^/^f{M-). Note that 
these operators are isometrics on L2(M'^). 

Definition 2.1. Let {i/jj : 1 < j < iV} C L2(M'^). The family 

* = {Ui^T^i/jj : 1 < i < TV, £ e Z, a e Z*^} 
is a wavelet tight frame of L2(M'^), if 

E UUi,T^^,)\' for all f e L,{R% (1) 

i<j<N,eei., 

The foundation for the construction of multiresolution wavelet basis or wavelet tight frame 
is a compactly supported real- valued function e L2(M'^) with the following properties. 
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(i) is refinable, i.e. there exists a finitely supported sequence p = {p{a))^^^d, G M, 
such that satisfies 

= m"^/2 ^ p{a)UMTa(i)ix). x E W^. (2) 
Taking the Fourier- Transform 

of both sides of (|2]) leads to its equivalent form 

^{M^uj) = P{e-''^)^{u), e-*"' = (e-*"'\ . . . , e''"^^), oo G (3) 
where the polynomial P E C[T] is given by 

For the characters x G G", we define the so-called isotypical components of P of 
weight X by 

Pje--)=m-i 5^ p(a)e"--, x G G'. (4) 

Q=X mod AfZ"* 

Clearly, P = E^^c ^x- 

(ii) One usually assumes that 0(0) = 1 by proper normalization. This assumption on and 
([3]) allows us to read all properties of from the polynomial P, since the refinement 
equation ([3]) then implies 

oo 

e=i 

The uniform convergence of this infinite product on compact sets is guaranteed by -P(l) = 
1. 

(iii) One of the approximation properties of is the requirement that the translates Tq,0, 
a G Z*^, form a partition of unity. Then 

P^(l)=m-\ xeG'. (5) 

The real- valued functions ipj, j = 1, . . . , N, are assumed to be of the form 

V;-(M^u;) = g,(e--)0H, 

where Qj G C[r] have real coefficients. These assumptions imply that ipj have compact support 
and, as in ([2]), are finite linear combination of UMTa<P- 
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We next describe the method called UEP (unitary extension principle) that allows us to 
determine the polynomials Qj, 1 < j < N, such that the family \E' as in Definition 12.11 is a 
wavelet tight frame of L2(M'^), see [HI [M]. To that purpose we need to introduce additional 
notation. Note that the group G acts on the coordinate ring C[T] by 

P ^ P"(e-^^) := P(e-''(^+")), a e G, w G 

and this operation commutes with the *— involution. 

Theorem 2.2. (UEP) Let the trigonometric polynomial P G C[T] satisfy P(l) = 1. If the 

trigonometric polynomials Qj G C[T], 1 < j < A^, satisfy the identities 

N 

S^,r - P"*(e-^")P"(e-^") = (e-^")g}(e-^'^), a, r G G, (6) 

i=i 

then the family is a wavelet tight frame o/L2(M°'). 

Certainly, the necessary condition for such Qj to exist is that the polynomial 1— P*P G R[T] 
is non-negative for all cj eW^. We give an equivalent formulation of the identities in Theorem 

Theorem 2.3. The identities are equivalent to 

U{e-'n*Uie-'n = Im, uje (7) 

with 

[/*=[P- ... Q-]^^^gM„x(7V+i)(C[T]). 

Moreover, this implies the "suh-QMF" condition 

1 - P"(e~*"')P"*(e-''^) > for all u G M^. (8) 

We note that ([H]) implies (iii). Moreover, the results in [52], see Theorem 12.41 imply that 
finding the sum of squares decomposition of / := 1 — 'Yla^G^'^^"* •> f ^ ^[^l^ is sufficient 
for the existence of the polynomials Qj in Theorem 12. 2[ The authors in [32] also provide a 
method for the construction of Qj from the sum of squares decomposition of the trigonometric 
polynomial /. 

Theorem 2.4. If the polynomial P G C[T] and the matrix-valued polynomial H G Mixs(C[T]) 
exist and satisfy 

^ p-(e-i<^)p-*(e-*^) + H{e~'^'^^'')H*{e-'^^^'^) = 1, (9) 

then there exist trigonometric polynomials Qi, . . . ,Qn G C[T] with real coefficients, N < m+s, 
that satisfy ([6]). 

Remark 2.5. Let Xj be the elements of G' = {xi, • • • , Xm} and Hj be the columns of if = 
{Hi,...,Hs) G Mixs(C[T]). The proof of Theorem 12.41 in [32] yields the explicit form of 
Qi, . . .,Qn, namely 

Q^.(e-*^) = m-i/2g-i^x. (i_^p(e-»-')p*^(e-*^))^ 1 < J < m, (10) 
Qm+Ae~n = P{e-nH,ie-'''"n, I < J < s. (11) 
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3 Equivalent formulations of UEP and moment prob- 
lems 

In this section we give several equivalent formulations of the Unitary Extension Principle stated 
in Theorem 12.21 that allow us to establish a connection between UEP and moment problems. 
Firstly, we reformulate UEP in terms of the isotypical components P^, Qj,x the polynomials 
P, Qj. By (cr, x) = e^'^'^, a E G and x ^ G', we denote the natural pairing between G and G'. 
Note that {a, x) is a root of unity of order dividing 

Theorem 3.1. Let P e C[T] satisfy P(l) = 1 and P^(l) = If e C[T], 1 <j <N, 

satisfy the identities 

N 



p;Px + J2QIxQj,x = ^-\ xeG', (12) 

N 

P;P, + J2 QlxQj,v = X,veG', Xt^V, 



then the family is a wavelet tight frame o/L2(M'^). 

Proof Note that P = ^ P^ and P^ = m'^ x)P'' imply 

xeG' aGG 

p* = J2p; and p''* = j:ip;r = J2{^,x)p;. 

xeG' xeG' xeG' 

Thus, 

x,»7'eG' jjeG' xeG' 

for rj' = X + V- Similarly for Qj. Therefore, the system in Theorem 12.21 is equivalent to 

E X) [Px^x+v + E ^L^i'X+r? ) = ^afiKo^ V ^ G' , a E G, 
xeG' V j=i ) 

For fixed ri G G', this is a system of m equations indexed by a G G in m unknowns P^P^+n + 

N 

E] Qj,xQi.x+'7' X e G'. The corresponding system matrix A = ({a,x))aeG,x<^G' is invertible 

and v4"^ = m'M*. Thus, ([I3]) is equivalent to ([12]). □ 

To establish the connection between UEP and moment problems, we rewrite the above 
result in matrix form. For this purpose we need to introduce additional notation. Let a set X 
be a subset of Z"' containing {a G Z'^ : p{a) ^ 0}. We also define the orthogonal projections 
Ex G rI-^I^I^I to be diagonal matrices with diagonal entries given by 

E [a a\ = [ ^' « = ^ ^^''^ aeX 
^ I 0, otherwise. 
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Theorem 3.2. Let a row vector p = m ^[p{a) : a G X], a column vector x{e = 
[e-^"-"^ : aeXf and P {e-''^) = p ■ x(e-*") satisfy P(l) = 1, P^{1) = m-\ If Qj {e-'^) = 
Qj ■ X (e~"^) with row vectors of length \X\, I < j < N, satisfy the identities 

(e^-) (^dtagip) -p^p-J2 ^J^^' j (*^'") =^ X, ^ e G', (13) 
then the family is a wavelet tight frame o/Z/2(M'^). 

Proof Note that the definition of implies that P^ {e^''^) = p- E^-x {e"'"^) and Qj^^ (e"*'^) = 
q • ■ E^-x (e"^'^), j = l,...,N. Due to P-^(l) = m"^ we get 



x^ (e^^) E^diag(p)E^x (6^*"^) = m"! 



and 

(e^'^) E^dmg{p)Er,x (e"^-) =0, x ^- 
Therefore, f|T3|) is equivalent to f|T2l) . □ 



Define the matrices 

AT 



i? = diag(p) - p^p and ^ = ^gjq^.. (14) 

i=i 

Then the task of constructing tight wavelet frames can be formulated as the following problem 
of semi-definite programming: for a given R find a matrix O G M'"^' ^ 1"^' such that 



(e^-) E^-0- E,x (e-^-) =0, x, ^7 G G\ (15) 



and i? + O is positive semi-definite. If such matrix O G rI^I^I-^I does not exist increase the set 
X. 

To determine the frame generators, set 5* = i? + O and find by standard factorization 
techniques from linear algebra. 

Note that the identities ([T^ are equivalent to the following linear constraints on the null- 
matrices O 

O^.^^e-'^'^-") = 0, x,r/GG', 

or, equivalent ly, 

^Oa,Q+^ = 0, rG{/3-a : a,/3GX}. 

Example 3.3. To illustrate the concept of null-matrices, we consider first a very prominent 
one- dimensional example of a Daubechies wavelet. Let 

P {e-'^) =px {e-"^) , p=i[l + v/3 3 + V3 3- V3 1- V3], 
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: a €{0,1,2,3}] . In this case M 



and X (e = [e * 
{0, 1}. Therefore, the orthogonal projections E-^ G 



2, G = {0,7r} and G' 



p4x4 



X G G', are given by 



= diag[l, 0,1,0] and = diag[0, 1, 0, 1]. 



By f[T4|) . we have 



R 



1 

64 



-2v^ 
2 



which is not positive semi-definite. Define 

-8^3 



= — 
64 








4v^ 


-2^3 




2 " 


2^3 


-6 




2^ 


-6 


12-2^3 


-6- 


M\/3 


2v^ 


-6 + 4^3 


4- 


-6^3 . 


8^3 





" 




-8^3 













8^3 -i 


^\/3 







-8^3 i 


^^3 . 





satisfying (fT5|) . Then S = R+0 is positive semi-definite, of rank one, and yields the well-known 
Daubechies wavelet, see [16] defined by 

(e"*") = l[^-V3 -3 + V3 3 + V3 -l-v^l - (e-^'^) . 
8 

Another two-dimensional example of one possible choice of an appropriate null-matrix 
satisfying (fT5|) is given in Example 14.41 

Another, very similar, way of working with null-matrices was pursued already in [0]. We 
first sketch the derivation of another equivalent formulation for UEP, see P, Theorem 2.4], 
keeping the notation as close as possible to the one in j6]. Multiplying the equations by 



m 



e ^'^^ and summing over a G G we get an equivalent system of equations 



N 

j=i 

This equivalent form of the UEP allows us to use semi-definite programming for tight 
wavelet frame construction. The null- matrices O are now part of the matrix S in the following 
result. 

Theorem 3.4. Let the sequences be given by 

J2 Cx(a)e~^"-" := - P (e"^'^) e'^'-'^pj^^^^^, x e G'. 

Then the following statements are equivalent: 

(i) There exist Qj G C[T], j = 1,...,N, satisfying UEP and whose coordinate degrees are 
determined by X. 
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(a) There exists a positive semi-definite matrix S = {sa,p)a,p<^x such that {—S,t = 0) solves 
the optimization problem 

minimize t 

subject to linear constraints 

Cx(«) = Sa+x+Mi3,x+Mp, X ^ G", a E {\ - fi : \, fi E 1} , 

and t-I-{-S)yO. 
For an application of the method described in Theorem 13.41 see [B]. 

4 Real algebraic geometry and wavelet frames 

In this section we use the result of Theorem 12. 4l that allows us to reduce the problem of existence 
of Qj in to the existence of an sos decomposition of a single nonnegative polynomial 

/ := 1 - ^P'^P'^* G M[T]. (16) 

crGG 

In subsection 14. H for dimension d = 2, we show that the polynomials Qi, . . . , Qn G C[T] as in 
Theorem 12.21 always exist. This result is based on recent progress in real algebraic geometry. In 
subsection 14. 2^ we give sufficient conditions for the existence of the Q/s in the multidimensional 
case. 

We start by deriving an equivalent formulation of the sub-QMF condition (|8]) that follows 
from part (a) of the following lemma. 

Lemma 4.1. Let P G C[T] with isotypical components P^, X ^ G' . 

(a) J2P-P-* = \G\-J2PxPx- 

<T&G xeG' 

(b) Given a decomposition 

r 

l-^P'^P'^* = ^HjH* (17) 

a&G j = l 

of length r with Hj G C[T]), then there also exists a decomposition 

r 

<7eG 3=1 xeG' 

of length mr with G-invariant matrix polynomials Hj^^ G C[T]. 

Proof, (a) follows from [H Lemma 2.3]. For (6) we observe that the left-hand side of (fT71) is 
G-invariant, i.e., 



o-eG ' ' j=i ueG 



Using the result in (a) we get that the left-hand side of (fT7|) is equal to 



I I j=l <7GG 3=1 xeG' 

For every x G G" choose a lift G Z'^ with respect to the natural surjection if — )■ G", a i— )■ Xo, 
given by 

Let ifj^^ := z~'^^Hj -^. Then, for all j and x, the polynomial Hj^^ is G-invariant and satisfies 



4.1 Scalar 2-dimensional case 

In this section we show that in the two-dimensional case {d = 2) the question of existence of 
a wavelet tight frame can be positively answered using the results from [35]. Thus, the main 
result of this section answers a long standing open question for the existence of tight wavelet 
frames. We also show that the result of Lemma 14. combined with Theorem 12. 4[ leads to an 
elegant method for the construction of tight frames. 

Theorem 4.2. Let d = 2, P E C[T] be a polynomial satisfying P{1) = 1 and X^^-gc < 1 
on = T(M). Then there exist N eN and polynomials Qi, . . . , G C[T] that satisfy 

N 

S^,r = P"*P^ + QTQh (y.reG. (19) 

Proof. The torus T is a non-singular affine algebraic surface over M, and T(M) is compact. The 
polynomial / in (ITB]) is in M[T] and is nonnegative on T(M) by assumption. By Corollary 3.4 
of there exist Hi, . . . ,Hr G C[T] satisfying / = X]j=i l-^iP- According to Lemma H?T] part 
(6), the polynomials Hj can be taken to be G-invariant. Thus, by Theorem 12. 4[ there exist 
polynomials Qi, ■ ■ ■ ,Qn satisfying ( 1T9|) . □ 

Lemma [4.11 part (a) sometimes yields an elegant method for the construction of Hj. We 
have 

So it suffices to find an sos decomposition for each of the polynomials — {Pxl"^, provided 
that they are all nonnegative. This nonnegativity assumption is satisfied, for example, for the 
special case when all coefficients p{a) of P are nonnegative. This is due to the simple fact that 
for nonnegative p{a) we get 

\Px\' < l^x(l)P = 

on T'^, for all x e G'. 

Example 4.3. As an illustration of the case of nonnegative coefficients we consider the three- 
directional piecewise linear box spline with the symbol 

Pizi,Z2) = lil + zi){l + Z2)il + ZiZ2), zj = e-'^^. 
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The SOS decomposition for the isotypical components yields 

/ = l-m^|P^p = isin2(a;i) + ^sin^(cu2) + ^sin^(a;i + a;2) 



Thus, in f lT7|) we have a decomposition with = 3. Since each of Hi, H2, H3 has only 
one isotypical component, we get a representation / = Hf + if| + if| with 3 G-invariant 
polynomials Hj. By Theorem 12.41 we get 7 frame generators. Note that the method in [32| 
Example 2.4] yields 6 generators, one of slightly larger support. The method in [HI Section 4] 
based on properties of the Kronecker product leads to 7 frame generators whose support is the 
same as the one of P. One can also employ the technique discussed in |2I] Section] and get 7 
frame generators. 

Another prominent example from bivariate subdivision is the so-called butterfly scheme. 
This example shows that the real advantage of treating the isotypical components of P sepa- 
rately is when one works with P of larger support. 

Example 4.4. The butterfly scheme describes an interpolatory subdivision method which is 
used in order to generate a smooth regular surface interpolating a given set of points [20] • The 
trigonometric polynomial P associated with the butterfly scheme is 

P{Zi, Z2) = ^ + ^ (^1 + ^2 + Z1Z2 + z{^ + Z2^ + Zi^Z2^^ 

+ ^ {zlz2 + Zizj + ZiZ2^ + Zi^Z2 + Z^"^ Z2^ + Z^^ Z2'^^ 

~64 (^^^^ ^^^^ ~^ ^^^^ ~'~ ^^^^ ~'~ ^^^2 ^ ~'~ ^^^"^ ^ 

-\-Zi ^Z2 + Zi Z2 + Z^ ^Z2 ^ -\- Z^ ^Z2 -\- Z^ Z2 ^ -\- Zi^ Z2^^ ■ 



Its first isotypical component is Po,o = |) which is the case for every interpolatory subdivision 
scheme. The other isotypical components, in terms of z^ = e~*'^*= {k = 1,2) are Pifi{zi, Z2) = 
\ cos(ci;i) + ^ cos(c<;i + 2^2) — ^ cos(3a;i + 2uj2) — ^ cos(a;i — 2^2), i.e., 

-^1,0(2:1,22) = 7Cos(c<;i) + ^sin^(wi)cos(wi + 2^2), 
4 o 

and Po,i(^i,2:2) = Pi,o(22,2i), Pi, 1(2:1, 2:2) = Pifi{ziZ2, Z2'^). Note that |P^| < \ for all x ^ G' . 
Elementary computation shows that 1 — 16 |Pi,o(-2:i, -22) | = 1— cos^(ci;i)— cos(a;i) sin^(ci;i) cos(ci;i-|- 
2^2) — I sin'^(ci;i) cos^(a;i + 2ci;2). Setting Uj := sin(ci;j), j = 1,2, v := sin(a;i +002), v' := 
sin(a;i — UJ2), w := sin(a;i + 2U2) and w' := sin(2c(;i + UJ2), we get 

1- 16|Pi,o(2i,22)r = ^ul(w^ + {ul + v^f + 2ul + 2v^y 

Therefore, 

i-^ip-|2 = l(uiui + uy + uy') + ^(uiw^ + uiw'^ + v%'^ 



16 
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This provides a decomposition 1 — Xlo-eG l-^^^P = Sj=i l-^il^ ^ ^ squares. As in 

the previous example, each Hj has only one nonzero isotypical component Hj -^_.. Thus, part 
(6) of Lemma 14.11 and by Theorem 12. 4[ there exists a tight frame with 13 generators. Namely, 
by Remark 1^751 we get 

11 1 

Ql{Zl, Z2) = - - 2^(^l' ^2), Q2{Zi, Z2) = -Zi - 2P{Zi, Z2)P^-^^o){Zi, Z2) 
Q?.{ZI,Z2) = Q2{z2,Zi), Qa{Zi,Z2) = Q2{ziZ2,Z2^) 
Qi+j{.Zi, Z2) = P{Zi, Z2)Hj^^., j = 1, . . . , 9, 

where Hj^^^ are the lifted isotypical components defined as in Lemma I^TTl Let X = {0, . . . , 7}^, 
P{zi,Z2) = p ■ x{zi,Z2) and Qj{zi,Z2) = q. ■ x{zi,Z2) with x{zi,Z2) = [z'^ : a E X]^. The 



corresponding null-matrix O E 



p 64x64 



defined in section [3] is given by 



x'^{z-^ ^ , Z2^) ■ O ■ x{zi,Z2) = X^{z^ \ -^2 



13 



^ gj^j - diag(p) + p^p 

.i=i 



X[Zi,Z2) 



Note that other factorizations of the positive semi-definite matrix diag(p) — p^p + O of rank 
13 lead to other possible tight frames with at least 13 frame generators. An advantage of using 
semi-definite programming techniques is that it can possibly yield Qj of smaller degree and 
reduce the rank of diag(p) — p^p + O. 

Using the technique of semi-definite programming the authors in [6] constructed numerically 
a tight frame for the butterfly scheme with 18 frame generators. The advantage of our present 
construction is that the frame generators are determined analytically. The disadvantage is 
that their support is approximately twice as large as that of the frame generators in [6J. 

The next example is one of the family of interpolatory -\/3— subdivision studied in 

[1 2 1 

associated dilation matrix is M = ^ ^ and m = 3. 
Example 4.5. The symbol of the scheme is given by 

^(^1,2:2) = P{0fl){Zx,Z2) + P(i,o) (21,^2) + P(0,l)(^l,^2) 

with isotypical components -P(o,o) 



The 



1 

3' 



^{Z2 + Z^^ + ZxZ2 ^) - ^(^1 '^Zl + zf + Z2 



P{0,l){Zl, Z2) 

and P(^i^o){zi, Z2) = P(o,i)(-22, -^i)- We have by Lemma I^TTl and due to |P(o,i)P 



|f(i,o)l 



creG 



(0,1)1 



thus it suffices to consider only 

m'^ - |Po,ir = 3"^-27"2(^51 + 16cos(u;i + W2) + 16cos(2wi-u;2) 

+ 16 cos((X'i — 2^2) + 2 cos(2(X'i + 2ijJ2) + 2 cos(2ci;i — 4c<;2) + 2 cos(4a;i — 2c<;2) 
— 8 cos(3(X'i) — 8 cos(3u;2) — 8 cos(3u;i — 3(X'2 j . 

Numerical tests show that this polynomial is nonnegative. 
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4.2 Scalar multivariate case 

In the general multivariate case > 2, in Theorem \A.7\ we provide a sufficient condition for 
the existence of a sums of squares decomposition of / in (ITB]) . This condition is based on the 
properties of the Hessian of / G M.[T] 

Hess(/)(e-) = (l^^/(e-)),,,.,H.,, 
where D'^ denotes the th partial derivative with respect to w G M'^. 

Theorem 4.6. Let V be a non-singular affine M.-variety for which V{M.) is compact, and let 
f G MlV] with f > on V(R). For every ^ G V(M.) with /(^) = 0, assume that the Hessian of 
f at ^ is strictly positive definite. Then f is a sum of squares in M.[V]. 

Proof. The hj^otheses imply that / has only finitely many zeros in Therefore the claim 

follows from [36], Corollary 2.17 and Example 3.18. □ 

Theorem 14.61 implies the following result. 

Theorem 4.7. Let P G C[T], and write f = I - J2aeG l^''!^- Assume that P(l) = 1 and 
f > on T(R) = T'^. If the Hessian of f is positive definite at any zero of f in T'^, there exist 
N eN and polynomials Qi, ■ ■ ■ ,Qn G C[T] satisfying ([6]). 



Proof. By Theorem 14.61 / is a sum of squares in M[T]. The claim follows then by Theorem 

1121 □ 

Due toP(l) = l, z = lis obviously a zero of /. We show next how to express the Hessian 
of / at 1 in terms of the gradient VP(1) and the Hessian of P at 1, if P additionally satisfies 
the so-called sum rules of order 2, or, equivalently, satisfies the zero conditions of order 2. We 
say that P G C[T] satisfies zero conditions of order k, if 

D'^P(e-'") = 0, /iGNjJ, \l2\<k, aeG\{0}, 

see [23 ES] for details. The assumption that P satisfies sum rules of order 2 together with 
P(l) = 1 are necessary for the continuity of the corresponding refinable function 0. 

Lemma 4.8. Let P G C[T] with real coefficients satisfy the sum rules of order 2 and P(l) = 1. 
Then the Hessian of f = 1 — J2aeG l-P'^P ^ equal to 

-2Hess(P)(l) - 2 VP(1)*VP(1). 

Proof. We expand the trigonometric polynomial P : T"' — )■ C in a neighborhood of 1 and get 

P(e-'^'^) = 1 + VP{l)u + ^w^Hess(P)(l)u; + 0(1^1^). 

Note that, since the coefficients of P are real, the row vector v = VP(1) is purely imaginary 
and Hess(P)(l) is real and symmetric. The sum rules of order 2 are equivalent to 

P"(l) = 0, VP"(1) = for all aeG\ {0}. 

Thus, we have P°"(e~"^) = for all a G G \ {0}. Simple computation yields 

|P(e-'''^)|2 = l + {v + v)uj + u^{Ress{P){l) + v*v)u + 0{\uf) 
= l + uj^{Ress{P){l)+v*v)uj + 0{\ujf). 

Thus, the claim follows. □ 
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Remark 4.9. Note that Hess(/) is a zero matrix, if P is a symbol of interpolatory subdivision 
schemes, i.e., 

xeG'\{0} 

and satisfies zero conditions of order at least 3. This property of Hess(/) follows directly from 
the equivalent formulation of zero conditions of order fc, see p]. The examples of P with such 
properties are for example the butterfly scheme in Example 14.41 and the three-dimensional 
interpolatory scheme in Example 14.121 

For simplicity of presentation, we start by applying the result of Theorem 14.71 to the 
2— dimensional polynomial / from Example 14. 3 [ This example also motivates the statements 
of Remark 14.111 

Example 4.10. The three-directional box-spline from Example 14.31 is defined by the trigono- 
metric polynomial 

P{e-n = e-'^-^^-^^ cos (^) cos (^) cos (^^^^j , G 

Note that 

Therefore, as the trigonometric polynomial P satisfies sum rules of order 2, we get 

Thus, the Hessian of / at 1 is positive definite. 

To determine other zeroes of /, by Lemma 14.11 part (a), we can use either one of the 
representations 

f{e-n = 1- E n 

(je{0,7r}2 0g{o,i}2\{o} 

= \Y. (l-cos^(^-x)). 

xe{o,i}2 

It follows that the zeros of / are all the points u G vrZ^ and, by periodicity of / with period 
vr in both coordinate directions, we get that 

Hess(/)(e-*") = Hess(/)(1), w G ttZ^ 

is positive definite at all zeros of /. 

Remark 4.11. {€) The result of [5], Theorem 2.4] implies the existence of tight frames for mul- 
tivariate box-splines. According to the notation in [iHl p. 127], the corresponding trigonometric 
polynomial is given by 

14 



.00 



where S = {^^^\ . . . ,^^'^^) E Z''^" is unimodular and has rank d. (Unimodularity means that 
all d X (i-submatrices have determinant 0, 1, or —1.) Moreover, S has the property that leaving 
out any column ^^^'^ does not reduce its rank. (This property guarantees continuity of the 
box-spline and that the corresponding polynomial P satisfies at least sum rules of order 2.) 
Then one can show that 

/(e— ) = l- V|p-(e--)|2>0, 



crGG 

and the only zeros occur at a; £ ttZ*^ and the Hessian of / at these zeros is positive definite. 
This yields an alternative proof for Theorem 2.4] in the case of box-splines. 
(ii) If the summands — |P^(e~"^)p are non-negative for all u e W^, then it can be easier 
to determine the zeros of / by determining the common zeros of all of these polynomials. 

Example 4.12. There was an attempt to define an interpolatory scheme for SD-subdivision 
with dilation matrix 2/3 in [14]. There are several inconsistencies in this paper and we give a 
correct description of the trigonometric polynomial P, the so-called subdivision mask. Note 
that the scheme we present is an extension of the 2-D butterfiy scheme to 3-D data in the 
following sense: if the data are constant along one of the coordinate directions (or along the 
main diagonal in M^), then the subdivision procedure keeps this property and is identical with 
the 2-D butterfiy scheme. 

We describe the trigonometric polynomial P associated with this 3-D scheme by defining 
its isotypical components. The isotypical components, in terms of Zk = e~*^*. A; = 1,2, are 
given by 

Po,0,o{^l,Z2,Z3) = 1/8, 

Pifi,o{zi,Z2,Z3) = -^coswi + ^(cos{uJi + 2^2) + cos(a;i + 2^3) + cos(wi + 2w2 + 2a;3) 
o 4 V 

A 



— ^cos(a;i — 2a;2) + cos(a;i — 2a;3) + cos(3a;i + 2a;2 + 2(^3) 

Po,l,o{zi,Z2,Z3) = Plfi,Q{z2,Zl,Z3), -^0,0,1 (^1 , ^2 , ^s) = -Pl,0,0 (^3 , ^1 , ^2) , 
-^1,1,1(^11^2,^3) = Plfi,o{ziZ2Z3,Z^'^,Z2^), 

Pi,i,o{zi,Z2,Z3) = - A^ cos(c<Ji +a;2) + A^cos(a;i -a;2) + cos(c<Ji +a;2 + 2a;3) 

^ cos(cJi -002 + 2(^3) + cos(a;i - UJ2 - 2^3) + 
cos(3cji + UJ2 + 2CJ3) + cos(a;i + 3uj2 + 2t^3)^ , 

Plfi,lizi,Z2,Z3) = Pl,lfiizi,Z3,Z2), Po,l,lizi, Z2, Z3) = Piflfi{z2, Z3, Zi), 

where A is the so-called tension parameter. 
The polynomial P also satisfies 

P{Zi, Z2, Z3) = hi + Zi){l + Z2){1 + Z3){1 + ZiZ2Z3)Q{zi, Z2, Z3) , Q{1) = 1, 
o 

which implies sum rules of order 2. 



15 



a) For A = 0, we have Q{zi, Z2, zs) = l/{ziZ2Zs). Hence, P is the scahng symbol of the 
trivariate box-spline with direction set (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 1, 1) whose support 
center is shifted to the origin. 

b) For < A < 1/16, the subdivision scheme converges and has a continuous limit function. 
The only zeros of the associated non- negative trigonometric polynomial / are at vrZ^, 
and the Hessian of / at these zeros is given by 

/ 1 - 16A i -8A I -8A \ 
Hess(/)(1) = Hess(/)(e-*"') = | - 8A 1 - 16A | - 8A , cue vrZ^. 

\ i - 8A i -8A 1 - 16A / 

The existence of the sos-decomposition of / is guaranteed by Theorem 14. 7t our explicit 
decomposition is computed as follows. 

bi) Denote u := cos(a;i+a;2), v := cos(a;i+a;3), w := cos(ci;2+Co'3), and u := sin(ci;i+Ci;2), 
V := sin(c<;i + u^), w := sin((X'2 + ^3). Elementary computations give 

and 

l-|Pi,i,oP = X\v'-w'f+[{^^ - Xv') + (1 - A.^^)) {lu'+X{v-uwr+X{w-uvr) 

which is SOS with 7 summands Hj, and each Hj has only one nonzero iso typical 
component. 

b2) The isotypical component -Pi,o,o is not bounded by 1/8; consider, for example, 
-Pi,o,o(e~*'^) at the point co = (— f , ~t)" obtain by elementary compu- 

tations 

-Pi, 0,0 = o coscJi + — Asino;!, A := sin2(ci;i + U2 + oJs) — sin2a;2 — sin2a;3, 
' ' 8 2 

and 

1 

16 
where 



|-Pl,0,or ~ |-Po,l,or ~ |-Po,0,lP — |-Pl,l,lP — -£^1,0,0 + -£^0,1,0 + -£^0,0,1 + -£^1,1,15 



„ 3A .4 /r> • A \2 1 — 16A .2/1 \ 

-c^ioo = T7:sm ui + —{Zsmui - AcosoJi) H — — sm Uiil + XA ); 

16 64 64 

the other -Ei.j,fc are given by the same coordinate transformations as Pij^k- Hence, 
we obtain an sos with 12 summands Hj, each with only one nonzero isotypical 
component. 

Thus, with the result of Theorem 12.41 we have explicitly constructed a tight frame with 
41 generators for the trivariate interpolatory subdivision scheme with tension parameter 
< A < 1/16. 
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c) For A = 1/16, the sum rules of order 4 are satisfied. In tliis particular case, the scheme 
is and the Hessian of / at 1 is the zero-matrix, thus the result of Theorem 14.71 is 
not applicable. Nevertheless, the sos decomposition of 1 — ^ \P"\'^ in b), with further 
simplifications for A = 1/16, gives a tight frame with 31 generators for the trivariate 
interpolatory subdivision scheme. Further properties of this scheme will be studied in 

The question may arise if there exists a trigonometric polynomial P that satisfies sum rules 
of a certain order and also satisfies the sub-QMF condition Y1icj(^g I-P'^P ^ 1; but there exists 
no UEP tight frame as in Theorem 12. 2[ In other words, can we find P as above such that the 
non- negative trigonometric polynomial 1 — |Pp is not an sos of trigonometric polynomials? 
Recall that such polynomials cannot exist in the 1-D and 2-D case. 

Theorem 4.13. There exists P G C[T] on the 3-D torus that satisfies sum rules of order 6 
and the sub-QMF condition, such that 1 — |Pp is not sos in M[T]. 

The proof is constructive. The following example defines a family of trigonometric polyno- 
mials with the properties in Theorem 14.131 We make use of the following local-global result: if 
the Taylor expansion of / G M[T] at one of its roots has, in local coordinates, a homogeneous 
part of lowest degree which is not sos of real algebraic polynomials, then / is not sos in ]R[T]. 

Example 4.14. Denote zj = 6"^'^^ j = 1,2,3. We let 

P{z) = (1 - cM{z))A{z), zeT, < c < -, 

3 

where 

M{z) = ytvl + ylvl - yl - 'iylylyl G M[T], y, = sinu;,. 

In the local coordinates (2/1,1/2, Z/s) at 2; = 1, M is the well-known Motzkin polynomial in 
K [2/1, 1/2, 2/3]; i-e. M is not sos in ]R[?/i, ?/2, 2/3]- Moreover, A G ]R[T] is chosen such that 

DM(1) = (5o,„, L'"A((t) = 0, < |a| < 8, a G G\ {1}, (20) 

and < 1. Such A can be, for example, any scaling symbol of a 3-D orthonormal 

wavelet with 8 vanishing moments; in particular, the tensor product Daubechies symbol A{z) = 
ms{zi)mg{z2)ms{z3) with mg in [16] satisfies conditions ( 120|) and = 1. The properties 

of M and A imply that 

1. P satisfies the sub-QMF condition, since M is G-invariant and < 1 — cM < 1, 

2. P satisfies sum rules of order at least 6, 

3. the Taylor expansion of 1 — at z = 1, in local coordinates (?/i, ?/2, 2/3), has 2cM as its 
homogeneous part of lowest degree. 

Therefore, 1 — P^ is not sos of trigonometric polynomials in M[T]. 

Next, we give an example of a trigonometric polynomial P satisfying P(l) = 1 and sum 
rules of order at least 2, but its corresponding polynomial / is negative for some w G M^. 
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Example 4.15. Consider 



P{Z1,Z2,Z3) = 6Z1Z2Z3 



3 /I , \ 3 /I , \ 3 



5 fl + zA fl + Z2y fl + ZsV fl + ZiZ2Zs\ 
5 fl + ZiV fl+Z2\ fl + ZsV fl + ZiZ2Z3Y 

[-^) [-^) [-^) [-^) - 

5 f^ + ZlV fi + Z2Y fi + ZsV fl + ZiZ2Zs\ 

The associated refinable function is in C(M"^) as the corresponding subdivision scheme is uni- 
formly convergent, but P does not satisfy the sub-QMF condition, as 

l-J2\P''ie-n\'<0 for a;=Q,0,0). 

o-eG 
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